Bidirectional cooperative motion of myosin-II motors on actin tracks with randomly 

alternating polarities 
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The cooperative action of many molecular motors is essential for dynamic processes such as cell 
motility and mitosis. This action can be studied by using motility assays in which the motion of 
cytoskeletal filaments over a surface coated with motor proteins is tracked. In previous studies of 
actin-myosin II systems, fast directional motion was observed, reflecting the tendency of myosin 
II motors to propagate unidirectionally along actin filaments. Here, we present a motility assay 
with actin bundles consisting of short filamentous segments with randomly alternating polarities. 
These actin tracks exhibit bidirectional motion with macroscopically large time intervals (of the 
order of several seconds) between direction reversals. Analysis of this bidirectional motion reveals 
that the characteristic reversal time, r rcv , does not depend on the size of the moving bundle or on 
the number of motors, N . This observation contradicts previous theoretical calculations based on 
a two-state ratchet model (Badoual et al., Proc. Natl. Acad. Sci. USA, 2002, 99, 6696), predicting 
an exponential increase of r rcv with N. We present a modified version of this model that takes into 
account the elastic energy due to the stretching of the actin track by the myosin II motors. The 
new model yields a very good quantitative agreement with the experimental results. 

PACS numbers: 



1. INTRODUCTION 

Cells utilize biological motors for active transport of 
cargo along their respective filaments to specific desti- 
nations [l||. Various types of motor proteins have dif- 
ferent preferred directions of motion. Most kinesins and 
myosins, for instance, move towards the plus end of mi- 
crotubules (MTs) and actin filaments, respectively Q. 
Others, such as Ned and myosin VI, move towards the 
minus end H,0| ■ While some processes, such as the trans- 
port of cargoes is achieved mainly by the action of indi- 
vidual motors, other processes, such as cell motility and 
mitosis, require the cooperative work of many motors. 
Muscle contraction, for instance, involves the simulta- 
neous action of hundreds of myosin II motors pulling 
on attached actin filaments and causing them to slide 
against each other Similarly, groups of myosin II 
motors are responsible for the contraction of the con- 
tractile ring during cytokinesis [f|. In certain biological 
systems, cooperative behavior of molecular motors pro- 
duces oscillatory motion. In some insects, for instance, 
autonomous oscillations are generated within the flight 
muscle 0. Spontaneous oscillations have also been ob- 
served in single myofibrils in vitro Q. Finally, dynein 
motors could be responsible for the oscillatory motion of 
axonemal cilia and flagella [E [l(| • 

The directionality of individual motors stems from in- 
teractions between different parts of the motor and from 
interactions between the motor and the track filament 
[TTI . [l2| . The direction of motion of a large collection of 
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motors may also be influenced by their cooperative mode 
of action. Specifically, in several recent experiments the 
ability of motors to cooperatively induce bidirectional 
motion has been demonstrated. These in vitro experi- 
ments were performed by using motility assays in which 
a filament glides over a dense bed of motors. In one such 
experiment, unidirectional motion of actin filaments due 
to the action of myosin II motors was transformed into 
bidirectional motion by the application of an external 
stalling electric field [HI]. Under such conditions the ex- 
ternal forces acting on the actin filament nearly balance 
the forces generated by the motors. Electric field was also 
used to bias the direction of motion in kinesin-MT sys- 
tems flij ]. In another experiment, bidirectional motion 
of MTs was observed when subjected to the action of an 
ensemble of NK11 motors [l5|. These motors are a mu- 
tant form of the kinesin related Ned, which individually 
exhibit random motion with no preferred directionality 
(TH - More recently, the motion of MTs on a bed of a 
mixed population of plus-end (kinesin-5 KLP61F) and 
minus-end (Ned) driven motors was shown to exhibit 
dynamics whose directionality depends on the ratio of 
the two motor species, including bidirectional movement 
over a narrow range of relative concentrations around the 
"balance point" [1 61 ] . Similarly, bidirectional transport of 
microspheres coated with kinesin (plus-end directed) and 
dynein (minus-end directed) on MTs was also reported 

m 

Several aspects of cooperativity in molecular motor 
systems have been addressed usin g d ifferent theoretical 
models pj, M, ES HE 11 11 H, M SS HI . One fea- 
ture which has not been treated in these studies is the 
dependence of the motion on the number of acting mo- 
tors. A notable exception is the work of Badoual et al. 
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[20j], where a two-state ratchet model has been used to 
examine the bidirectional motion observed in the NK11- 
MT motility assay described in ref. 15]. The model of 
Badoual et al. [2(| demonstrated the ability of a large 
group of motors working cooperatively to induce bidi- 
rectional motion, even when individually the motors do 
not show preferential directionality. (The model also pre- 
dicts that directional motors can also induce bidirectional 
movement, if the filaments are close to stalling condi- 
tions in the presence of an external load.) According 
to this model, the characteristic time in which the fila- 
ment undergoes direction reversal ("reversal time"), r rcv , 
increases exponentially with the number of motors, N. 
Thus, the reversal time diverges in the "thermodynamic 
limit" N — > oo, and the motion persists in the direction 
chosen at random at the initial time. 

In this work, we present an in vitro motility assay in 
which myosin II motors drive the motion of globally a- 
polar actin bundles. These a-polar bundles are generated 
from severed (polar) actin filaments whose fragments are 
randomly recombined. When subjected to the action of 
a bed of myosin II motors, these a-polar bundles exhibit 
bidirectional motion with characteristic reversal times 
that are in the range of T rev ^3 — 10 sec. The reversal 
times of the dynamics show no apparent correlation with 
the size of the gliding bundles, or equivalently, with the 
number of motors N interacting with the track (which, 
because of the homogeneous spreading of the motors on 
the bed, is expected to be proportional to the size of the 
moving bundle). This observation is clearly in disagree- 
ment with the strong exponential dependence of r rev on 
N, predicted by Badoual et aZ. [20j . 

Here, we propose a modified version of this model 
that explains the experimentally observed independence 
of r rcv on N. We argue that the origin of this behavior 
can be attributed to the tension developed in the actin 
track due to the action of the attached myosin II motors. 
An increase in the number of attached motors leads to an 
increase in the mechanical load which, in turn, leads to an 
increase in the detachment rate of the motors, as alr eady 
suggested in models of muscle contraction (23.l23.l25ll2q|. 
Unlike most previous studies where the myosin conforma- 
tional energy was calculated, in this work we consider the 
elastic energy stored in the actin track and demonstrate 
that the detachment rate increases exponentially with N. 
This unexpectedly strong effect (which is another, indi- 
rect, manifestation of cooperativity between the motors) 
suppresses the exponential growth of r rev with AT. 



2. MATERIALS AND METHODS 
A. Protein purification 

Actin was purified from rabbit skeletal muscle acetone 
powder p7j . Purification of myosin II skeletal muscle is 
done according to standard protocols [28j . Actin labeled 
on Cys374 with Oregon Green (OG) purchased from In- 



vitrogen. 

B. NEM myosin II 

N-ethylmaleimide (Sigma, Co.) inactivated myosin II 
was prepared according to standard protocol of Khun and 
Pollard HI]. 

C. Optical Microscopy 

Actin assembly was monitored for 30 minutes by fluo- 
rescence with an Olympus IX- 71 microscope. The labeled 
actin fraction was 1/10 and the temperature at which the 
experiments were conducted was 23° C. Time-laps images 
were acquired using a DV-887 EMCCD camera (Andor 
Co., England). 

D. Motility assay 

Protocol for this assay was adopted from Kuhn et al. 
[29l ] . The assay includes two essential steps: (a) immo- 
bilization of actin filaments on a bed of NEM myosin II 
inactivated motors, and (b) addition of active myosin II 
motors at a defined concentration. For that purpose, 
7.5-8.5/il of Q.2/iM NEM myosin II is introduced into 
a flow chamber (26mm x 2mm glass surface area) for 1 
minute of incubation followed by wash of the flow cham- 
ber with BSA solution to passivate the surface. Following 
this, actin filaments were grown on the surface (3/iM 10% 
O.G. labeled). Finally, the cell was supplemented with 
8/A of 0.6/xM myosin II motors (in 2X myosin solution 
containing: 3.3mM MgCl 2 . 2mM EGTA, 20mM HEPES 
pH=7.6, 1% MethylCellulose, 3.34mM Mg-ATP, 400mM 
DTT, 17.6mM Dabco), supplemented with 0.133M KC1, 
5/jM Vitamin D, and an ATP regenerating system con- 
taining O.lmg/ml Creatine Kinase (CK) and ImM Crea- 
tine Phosphate (CP). At the KC1 concentrations used in 
this assay, the myosin II motors are assembled in small 
motor aggregates (~16 myosin II units/aggregate) also 
known as mini-filaments [301 ]. Fluorescent images were 
taken every 2 seconds for 30 minutes. 

E. Data Analysis 

The position of fluorescent bundles was determined 
as the intensity center of mass using METAMORPH 
(Molecular Devices) software. The position was analyzed 
using a custom MATLAB (The MathWorks, Inc.) pro- 
gram. The data was corrected for stage drift. We first 
measured the fluctuations of the positions of the bun- 
dles in the absence of ATP (i.e., when the motors are 
not active). Under such conditions, the positions of the 
bundles measured every 1 sec exhibit a Gaussian distri- 
bution with zero mean and standard deviation A ~ 200 
nm (Supplementary Information ? Figure SI). Bidirec- 
tional motion (with ATP) was evaluated based on snap- 
shot taken every 2 sec. We set A as the experimental 
uncertainty since position changes that are smaller then 



A cannot be unambiguously identified as a change in the 
direction of motion caused by the action of the motors. 

F. Estimation of number of interacting motors 

In order to evaluate the number of acting motors, the 
dimensions of the bundles and the motors surface concen- 
tration, C m , must be determined. We estimate C m by 
assuming that all the motors that were introduced into 
the flow chamber adhere to the top and bottom glass 
surfaces of the flow cell (total surface, 104 mm 2 ). This 
gives C m ~ 27800 (/j,m)~ 2 , which corresponds to densely 
packed motor beds (typical distance of a few nanome- 
ters between motor heads). At such high densities, inho- 
mogeneities associated with the assembly of motors into 
mini-filaments can be ignored. 

The length of a bundle, L, was measured using META- 
MORPH software. The width of a bundle was estimated 
by dividing its fluorescence intensity by the intensity of 
single actin filaments, which gives an estimate for the 
number of filaments, Nf, composing the bundle. Assum- 
ing that the shape of the bundle is cylindrical, its radius 
can be estimated as R = y/~Nf ■ r, where r ~ 3.75 nm is 
the actin filament radius. The motors can only interact 
with the part of the bundle that faces the myosin bed 
covered surface. Assuming that this part corresponds to 
roughly a quarter of the surface of the bundle, we find 
that the area that comes into contact with the motors 
A ~ (%/2)RL. The number of interacting motors is, 
thus, given by N — C m A. Using this approximation and 
the measured surface concentration and bundle dimen- 
sions, we estimate (see Fig. [5]) that N ~ 1000 — 5000 for 
most of the bundles studied in this work. 
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FIG. 1: (A) Schematic diagram of the system before addition 
of active motors. The surface of a microscope slide was sat- 
urated by BSA (blue balls) and NEM myosin II (long, two- 
headed, brown objects). Actin filaments/bundles (thin yel- 
low line) are attached to NEM myosin heads above the sur- 
face. (B,C) Images of the system before the addition of active 
myosin II minifilaments. (B) shows the thick actin bundles 
formed at a high concentration of MgCl2(1.67 mM), while in 
(C) the thin bundles/filaments formed at a low MgCb con- 
centration (0.5 mM) are shown. Bar size is 5 /im. 



were identified by analyzing the position of this motor. 
The distribution of reversal times, t, follows an expo- 
nential distribution: pit) — (l/r rov ) exp(— t/r lev ), from 
which t 1cv was extracted. The error bars in Fig. [7p 
represent one standard deviation of the distribution of 
reversal times measured for different realizations of glob- 
ally a-polar tracks of similar size. For each value of N, 
the number of simulated realizations is 40. 



RESULTS 



G. Computer simulations 

A detailed discussion on the computational model is 
found in the section [H below. The model is based on 
the model presented in ref. > where N rigidly coupled 
equidistant motors interact with a one-dimensional peri- 
odic potential representing the actin track. The spacing 
between the motors q is larger than and incommensu- 
rate with the periodicity of the potential, I. The track 
consists of M ~ (q/l)N periodic units, which are repli- 
cated periodically. In each unit of the track, a force of 
magnitude / ran and random directionality is introduced 
which defines the local polarity of the track. Globally a- 
polar tracks were generated by setting the total random 
force to zero (i.e., choosing an equal number of periodic 
units in which the random forces point to the right and 
left). The motion of the motors on these tracks was cal- 
culated by numerically integrating the equations of mo- 
tions [dx = (F tot /X)dt] (see Eq. Q] and following text) 
with time step dt = 0.05 msec. The position of one of 
the motors along the track was recorded every 0.25 sec, 
and changes in the direction of the motion of the motors 



A diagrammatic representation of the system is dis- 
played in Fig. [TJ The protocol, based on that of Kuhn 
et al. [29|, is described in detail above. In brief, the 
surface of a microscope slide was saturated with NEM- 
inactivated myosin II motors (drawn as long, two headed, 
brown objects at the bottom of Fig. [TJ\) and passivated 
by BSA (blue balls in Fig. [TK). Subsequently, actin fila- 
ments/bundles (thin yellow line, Fig.[T]A_) were grown and 
held firmly on the underside of the NEM-myosin II bed. 
Fig. QJ3 shows a characteristic fluorescent microscope im- 
age of the system which, at this stage, consisted of a 
large number of long actin bundles. The bundles were 
formed due to the presence of free Mg 2+ ions (concen- 
tration 1.67 mM), which induced attractive electrostatic 
interactions between the actin filaments (3l| . Unlike bun- 
dles formed by certain actin-binding proteins, filaments 
formed by condensation in the presence of multivalent 
cations are randomly arran ged within the bundles with- 
out any specific polarity [32I l33| . At lower concentrations 
of Mg 2+ (0.5?mM), both thinner bundles and single fila- 
ments were observed (Fig. HP). 

After the initial step, myosin II minifilaments (multi- 
headed brown objects, Fig. [2JA.) were added to the cell 
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FIG. 3: Sequence of snapshots showing actin bundles moving 
directionally on a bed of active myosin II motors, at a low 
motor concentration ([myosin II motor]= 0.3 /iM). As a visual 
guide, we have marked one such actin bundle in orange. This 
bundle moved directionally across the image plane from left 
to right at times (A) 0:00, (B) 1:40 min, and (C) 7:20 min 
(typically at a velocity of 10 ^m/min). Bar size is 5 jj,m. 




FIG. 2: (A) Schematic diagram of the system after addi- 
tion of active myosin II motors. After the initial step (see 
Fig. [TJ , myosin II minifilaments (multi-headed brown objects) 
were added to the cell sample. The motors that landed on 
the BSA surface created a homogeneous bed of immobile, 
yet active, motors. Other motors landed on the actin fila- 
ments/bundles (long yellow line) present on the surface. The 
myosin II minifilaments started to move along the actin fil- 
aments/bundles. During their motion, the motors exerted 
forces on the actin filaments, which caused severing of small 
actin fragments (short yellow lines) . The ruptured actin frag- 
ments could move rapidly on the bed of active myosin II 
minifilaments and fuse with other bundles. One fusion event 
is demonstrated in the sequence of snapshots (B-E). Here, we 
show (B) two bundles moving oppositely to each other, getting 
closer (C) and then fusing (D-E) to create one larger object. 
Time is given in minutes, bar size is 5 /im. (F) The bundles 
continue to grow in size through multiple fusion processes, un- 
til eventually a large, highly a-polar bundle is formed (thick 
yellow tube - the inset illustrates the internal structure of 
such a bundle, consisting of individual actin filaments with 
randomly orientated polarities). 



sample. The motors that landed on the BSA surface 
created a homogeneous bed of immobile, yet active, mo- 
tors. Other motors landed on the actin filaments/bundles 
present on the surface. These motors started to move 
along the actin tracks, thereby exerting forces on the 
actin filaments, which led to the severing of small actin 
fragments [U Hfl (Fig. [5]A.) . The ruptured actin frag- 
ments were then free to move rapidly on the bed of ac- 
tive myosin II motors. When gliding bundle fragments 



came into close proximity to other bundles, they could 
fuse, creating new, a-polar, bundles (Fig. See also 
Supplementary Information - movie 1, showing a small 
actin piece severed from a filament, moving rapidly and 
fusing with a distant existing bundle). These newly cre- 
ated bundles could further fuse with each other to form 
even larger objects (see Fig.lU 7 and the sequence of snap- 
shots in Figs. [2J3-E depicting one such event of fusion 
of bundles). The rate of fusion events decreased with 
time and, after several minutes, the system relaxed into 
its final configuration, shown schematically in Fig. HJ 7 ). 
Notice that the severing and rearrangement of the origi- 
nally formed actin filaments/bundles (Fig.[Tj3) led to the 
formation of much shorter bundles (Figs. 03-E). More- 
over, the random nature of the multiple fusion processes 
involved in the generation of these shorter bundles en- 
sured that the final actin tracks were highly a-polar. In- 
deed, the motion of most of the bundles shown in Supple- 
mentary Information - movie 2 was bidirectional ( "back 
and forth" motion), and only those bundles undergoing 
rare fusion events exhibited unidirectional motion. It is 
important to emphasize that bidirectional motion was 
observed only above a certain concentration of added 
myosin II motors (0.6 /jM) and only in the presence 
of ATP. At lower concentrations of motors (0.3 fjM), 
the motion of actin bundles was directional (see Fig. [3] 
and Supplementary Information - movie 3. Note that in 
movie 3, the motion takes place both along pre-existing 
actin tracks, as well as on the BSA bed, both which are 
covered by active myosin II mini-filaments. The motion 
between these two areas is continuous, demonstrating 
that the whole surface is covered uniformly with motors.) 
We, therefore, conclude that the bidirectional movement 
originates from the action of the active myosin II motors 
which (i) severed actin pieces, (ii) transported the severed 
fragments, which fused into actin tracks with randomly 
alternating polarities, and (iii) moved these a-polar actin 
tracks bidirectionally (see Supplementary Information - 
movie 2). 

Fluorescence microscopy was used to follow the bidi- 
rectional motion of the actin tracks. Fig. [3J\ shows the 
position of center of mass of one bundle (three snapshots 
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FIG. 4: Position of a bundle over a time interval of 800 sec. 
The time interval between the consecutive data points is 2 sec. 
(B-D) Pseudo-color images of the actin bundle. The yellow 
arrows indicate the instantaneous direction of motion of the 
bundle. Bar size is 5 /jm. 
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are shown in Figs.[3j3-D) during a period of more than 10 
minutes of the experiment. The dynamics of this bun- 
dle are representative of the motion of the other actin 
bundles. Specifically, the one-dimensional motion of the 
bundle does not persist in the initial direction, but rather 
exhibits frequent direction changes. Measurements of the 
position of the center of mass of the bundle were taken 
at time intervals of At = 2 sec, and the mean velocity in 
each such period of motion was evaluated by v = Ax/ At, 
where Ax is the displacement of the center of mass (see 
section [5]). Fig. [SJA. shows the velocity histogram of the 
bundle shown in Fig. 3J The velocity histogram is bi- 
modal indicating bidirectional motion. The speed of the 
bundle varies between |u| = 1 — 2 /im/min, which is 2 
orders of magnitude lower than the velocities measured 
in gliding assays of polar actin filaments on myosin II 
motors [36[ . The fact that the typical speed of the bidi- 
rectional motion is considerably smaller than those of 
directionally-moving polar actin filaments can be par- 
tially attributed to the action of individual motors work- 
ing against each other in opposite directions. The bidi- 
rectional movement consists of segments of directional 
motion which typically last between 2 to 10 time inter- 
vals of A< = 2. The statistics of direction changes is 
summarized in Fig. [5)3 which shows a histogram of the 
number of events of directional movement of duration t. 
The characteristic reversal time, r rov , can be extracted 
from the histogram by a fit to an exponential distribu- 
tion: p{t) = (1/T rcv ) exp(— t/r ICV ). This form (which, as 
exemplified in Fig. [SJ3, fits the data well) is expected if 
the probability per unit time to "turn" in the opposite 
direction is independent of the time since the beginning 



FIG. 5: Velocity histogram of the bundle whose motion is 
shown in Fig. [4] (based on 900 sampled data points ? Figure 
U shows only 400 of those points), exhibiting a clear bimodal 
distribution. (B) Distribution of the reversal time for the 
same bundle. The distribution is fitted by a single exponential 
decay function with a characteristic reversal time: r rcv ~ 3 
sec. 



of motion in a given direction. 

Although Figs. [4] and [5] summarize the results corre- 
sponding to the movement of a single actin bundle, these 
results are representative of several tens of bundles whose 
motion we followed in several repeated experiments. In 
these experiments we observed that essentially all a-polar 
bundles exhibited bidirectional motion. For the sake of 
our quantitative analysis, we picked a smaller group of 
19 bundles (about 25% of all bundels) for which both 
the reversal time and the number of acting motors N 
(which is proportional to the surface area, see section [2]) 
could be determined with sufficiently high precision. The 
results corresponding to the motion of this sample of rep- 
resentative bundles are plotted in Fig. [Sj The choice of 
which bundles to include in Fig. [5] is based on the follow- 
ing practical reasons: The fluoresces intensity of small 
bundles is too low (compared with the background) and, 
thus, the accurate position and dimensions are hard to 
determine. Very large bundles arc practically immobile 
and their motion is smaller than the experimental spatial 
resolution (see Material and Methods - section[2]). Thus, 
the data in Figure [B] includes only bundles of "interme- 
diate" size. One can see in Fig. [5] that while N varies 
over half an order of magnitude, the corresponding T rev 
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FIG. 6: Characteristic reversal time, r rcv , of 19 different bun- 
dles as a function of the number of working motors TV. The 
reversal time for each bundle is obtained by an exponential 
fit. 

are similar to each other (3 < r rcv < 10 sec) and show no 
apparent correlation with N. 

4. DISCUSSION 

The bidirectional motion of motors was previously seen 
in systems consisting of motors that lack specific direc- 
tionality lla. 13711 , mixtures of motors working in opposite 
directions [16|] , or under the action of external forces close 
to stalling conditions (forces acting on the filament that 
nearly balance the forces generated by the motors) [l3[. 
One commonly used model for the dynamics of molecu- 
lar motors in the biophysics literature is the Brownian 
ratchet mechanism 38] . Within this modeling approach, 
the motion of individual motor proteins is studied by 
considering the motion of a particle in a periodic, locally 
asymmetric, potential. It follows from the second law of 
thermodynamics that if the system is coupled to a ther- 
mal bath, the particle subjected to the periodic potential 
will not exhibit large scale directed motion. Directed mo- 
tion is possible only if the system is (i) locally asymmet- 
ric, and (ii) driven out of equilibrium by an additional 
deterministic or stochastic perturbation. This perturba- 
tion is used in the model to represent the consumption 
of ATP chemical energy by the motors. Ratchet models 



are not molecular in nature but rather present a way to 
identify the minimal physical requirements for the mo- 
tion of motor proteins. However, by choosing properly 
the parameters of system, they may be employed to de- 
rive quantitative predictions for specific motor-filaments 
systems. Ratchet models have been extended for describ- 
ing and analyzing the collective motion of groups of mo- 
tors. The motion of several motors is influenced by the 
motor-motor interactions [39| and mechanical coupling 
(20| . The model proposed by Badoual at al. [2(| (and 
which, below, we present a slightly modified version of) 
demonstrates that mechanical coupling between the mo- 
tors is sufficient for the generation of highly cooperative 
bidirectional motion, even if the motors attach to/detach 
from the track in an uncorrelated fashion. 

A key prediction of the model in Ref. [2(| is the ex- 
ponential increase of the mean reversal time of the bidi- 
rectional motion, r rov , with N, the number of motors. 
This prediction is in contradiction with our experimental 
results (Fig. [|5]). Here, we show that this disagreement 
can be resolved by considering the stretching energy in- 
volved in the interactions between the actin track and 
the "walking" motors. Accounting for this effect elimi- 
nates the exponential dependence of r rcv on N. More- 
over, when values representing myosin II-actin systems 
were assigned to the parameters of the model, we found 
T rev ~ 1 — 12 sec, which in a very good quantitative agree- 
ment with the experimental data. 

Our model is illustrated schematically in Fig. [7jA: We 
consider the ID motion of a group of N point particles 
(representing the motors) connected to a rigid rod with 
equal spacing q. The actin track is represented by a peri- 
odic saw-tooth potential, U(x), with period I and height 
H. We choose q — (5n/12)l ~ 1.309Z, which satisfies the 
requirements of the model [20] for q to be larger than 
and incommensurate with the periodicity of the poten- 
tial. The locally preferred directionality of the myosin II 
motors along the actin track is introduced via an addi- 
tional force of size / ran exerted on the individual motors. 
In each unit of the periodic potential, this force randomly 
points to the right or to the left (the total sum of these 
forces vanishes), which mimics the random, overall a- 
polar, nature of the actin bundles in our experiments. 

The instantaneous force between the track and the mo- 
tors is given by the sum of all the forces acting on the 
individual motors: 



N 

*y ^ ymotor 



N 

£ 

i=l 



dU ( Xl + (i-l)q) 
dx 



+ /ran (^1 + (i - 1) q) 



Ci(t), 



(1) 



where Xi = x\ + (i — l)q is the coordinate of the i-th motor. The two terms in the square brackets represent 
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FIG. 7: (A) N point particles (representing the motors) are connected to a rigid rod with equal spacing q. The motors interact 
with the actin track via a periodic, symmetric, saw-tooth potential with period I and height H. In each periodic unit, there is 
a random force of size /ran, pointing either to the right or to the left (red arrows). The motors are subject to these forces only 
if connected to the track. The detachment rate uii is localized in the shaded area of length 2a < I, while the attachment rate 
LU2 is located outside of this region. The off rate 1^3 is permitted only outside the gray shaded area. (B) A set of n + 1 point 
particles connected via TV identical springs with spring constant k. Each particle is subjected to a random force ft that takes 
three possible values: — /, 0, +/ (blue arrows). The mean force acting on the particles is given by / = /cm/(t> + 1), where 
/cm = X^flt 1 /* i s the total force acting on the center of mass and causing the movement of the system. The force Fi stretching 
the i-th spring (red arrow) is equal to the sum of the excess forces, fi — f, acting on the particles located to the right of the i-th 
spring, Fi — ($3}=i fj) ~ Because a similar expression can be written taking into account the forces acting on the particles 
located to the left of the spring, one can readily show that J3 i=1 Fi = 0. (C) The mean reversal time, r rov , as a function of the 
number of motors N, computed for different realizations with a = 0.0018 (solid circles and a dashed line as a guide to the eye) 
and a = (open circles). In the latter case, the effect of actin stretching is neglected, and r rGV grows exponentially with N 
(dash-dotted line), in agreement with ref. [201 ]. The half-filled circles denote the experimental results, also presented in Fig. [6] 
(D) The fraction of attached motors, N c /N, as a function of the total number of motors, for a — 0.0018. The values of N c /N 
for a — are indistinguishable. 



the forces due to the symmetric saw-tooth potential and 
the additional random local forces acting in each periodic 
unit. The latter are denoted by red arrows in Fig. [7^\. 
The function Cj(t) takes two possible values, or 1, de- 
pending on whether the motor i is detached or attached 
to the track, respectively, at time t. The group veloc- 
ity of the motors (relative to the track) is determined 
by the equation of motion for overdamped dynamics: 
v(t) = i<t t(i)/A. The friction coefficient, A, depends 
mainly on motors attached to the track at a certain mo- 
ment and is therefore proportional to the number of con- 
nected motors, N c < N at time t: A — XoN c . 

To complete the dynamic equations of the model, we 
need to specify the transition rates between states (0 - 
detached; 1 - attached). The motors change their states 
independently of each other. We define an interval of size 
2a < I centered around the potential minima (the gray 
shaded area in Fig.JTJ^). If located in one of these regions, 



an attached motor may become detached (1 — > 0) with a 
probability per unit time u>\. Conversely, a detached mo- 
tor may attach to the track (0 — > 1) with transition rate 
UJ2 only if located outside this region of size 2a. However, 
we also allow another independent route for the detach- 
ments of motors, which may take place outside the gray 
shaded area in FigOV (i.e., around the potential max- 
ima) and is characterized by an off rate L03. The rates 
u>i, u>2, (see blue arrows in Fig. 7A), represent the 
probabilities per unit time of a motor to (i) detach after 
completing a unit step, (ii) attach to the track, or (hi) 
detach from the track without completing the step. 

Generally speaking, the rates of transitions between 
states depend on many biochemical parameters, most no- 
tably the types of motors and tracks, and the concentra- 
tion of chemical fuel (e.g., ATP). They may also be af- 
fected by the forces induced between the motors and the 
filament, which result in increase in the configurational 
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energy of the attached myosin motors [HI, [24], [2f| [26j 
and in the elastic energy stored in the S2 domains of the 
mini-filaments, as well as an increase in the stretching 
energy of the actin filament. The latter contribution can 
be introduced into the model via a modified detachment 
rate given by: 0J3 — uj® exp(— AE/UbT), where AE is 
the change in the elastic energy of the actin track due 
to the detachment of one motor head. The dependence 
of AE on the number of connected motors N c (out of a 
total number of motors, N) can be estimated in the fol- 
lowing manner: Consider a scries of N + 1 point particles 
connected by N identical springs (representing a series 
of sections of actin filaments) having a spring constant k 
(see Fig. 7B). Let us assume that random forces act on 
the particles and denote the force applied on the particle 
with index i (1 < i < N + 1) by fa. Assume that each 
of these forces can take three possible values: — / (repre- 
senting attached motors locally pulling the track to the 
left), +/ (attached motors pulling the track to the right), 
and (detached motors not applying force). Defining the 
"excess force" with respect to the mean force acting on 
the particles: // = /(- / [where / = fi/(N + 1)], 

one can show that the force stretching (or compressing) 
the t-th spring in the chain is given by the sum of excess 
forces acting on all the particles located on one side of 
the spring 

i N+l 
j=l j=i+l 

From Eq. [2] it can be easily verified that Yj%=i = 0- 
We thus conclude that the excess forces acting on the 
particles, /*, represent a series of random quantities with 
zero mean. Therefore, the size of Fi can be estimated 
by mapping the chain of springs into the problem of a 
ID random polymer ring [401 ]. where the elastic energy 
stored in the i-th spring, e, = Ff/2k, plays the role of the 
squared end-to-end distance between the i + 1 monomer 
and the origin. From this mapping we readily conclude 
that the energy of most of the springs (except for those 
located close to the ends of the chain) scales linearly with 
the number of attached motors: e ~ N c (f 2 /2k). The 
total elastic energy of the chain scales as 

E~Ne~NN e (f*/2k), (3) 

and when a motor detaches from the track (N c — > N c — 1), 

AE/k B T = —olN (4) 

where a is a dimensionless pref actor. 

We simulated the dynamics of an iV-motor system, 
choosing parameters corresponding to the myosin Il-actin 
system. The period of the potential / = 5 nm corre- 
sponds to the distance between binding sites along the 
actin track (4l|, E2, EH , and the amplitude of the sym- 
metric potential is set to H = 6k B T. Thus, the force 
generated by each motor head on the track is 2H/1 = 10 



pN (first term in square brackets in Eq. [T]) . The magni- 
tude of the random force that defines the local polarity 
of the track (second term) is given by / ran = 4.5 pN, so 
the total force acting on each motor head ranges between 
about 5 to 15 pN [5], El, El El • The interval around 
the potential minima from which motors can detach from 
the track with rate lo\ is chosen to be 2a = 3.8 nm. The 
transitions rates between attached and detached states 
are 1 = 0.5 ms and u^ 1 = 33 ms 0, E3, El, El ■ With 
this choice of parameters, we obtain a system with a low 
fraction of attached motors Nc/N ~ 0.1 (see Fig. 03). 
We also set the friction coefficient per attached motor to 
Ao = 1.25 • 10~ 4 kg/s, which yields the experimentally 
measured velocity v ~ 0.03 nm/ms ~ 2 /zm/min. The 
rate W3 expresses the probability of a single motor head 
to detach from the track without advancing to the next 
unit. The probability p of such an event is 1-2 orders of 
magnitude smaller than the complementary pro bability 
(1— p) to execute the step. We take p = 1/30 [49], which 
yields l^®)^ 1 ~ pv/l — 7500 ms. Finally, the exponent 
a appearing in Eq. 3] is evaluated by: 

a ~ (f/2kk B T) = (fl/2YAk B T), (5) 

where Y ~ 10 9 Pa is Young's modulus for actin and 
A ~ 35 nm 2 is the cross sectional area of an actin filament 
49]. For the model parameters: / ~ 10 pN, I = 5 nm, 
we find a ~ 0.0018. 

Fig. [7p shows the computationally measured rever- 
sal time r rev function of N for 800 < N < 3600. 
This range largely overlaps with the estimated range of 
number of motors in our experiments (see Fig. [SJ The 
experimental data points in this range of N are replot- 
ted in Fig. [7p and denoted in half-filled circles). For 
each N, the computational results represent the aver- 
age r rov computed for 40 different realizations of ran- 
dom, overall a-polar, tracks. The error bars represent 
the standard deviation of r rev between realizations, where 
for each realization r rcv is estimated by fitting the his- 
togram of turning times to an exponential decay func- 
tion (as in the experimental part - see Fig. 03). Two sets 
of computational data are shown in Fig. [7p: one cor- 
responding to a = 0.0018 (solid circles), and the other 
to a — (open circles), i.e., without considering the ef- 
fect of actin stretching, but when all the other system 
parameters mentioned above are kept unchanged. The 
latter case is qualitatively similar to the model presented 
in ref. [20], exhibiting a very strong exponential depen- 
dence of T rev on N (indicated by the straight dashed- 
dotted line in Fig. [7p). In contrast, the data correspond- 
ing to a — 0.0018 show much weaker variation in r rcv 
upon changing N. The mean reversal times computed for 
1200 < N < 2800 are found in the range 2 < r rcv < 12 
sec, which is in a very good quantitative agreement with 
the corresponding range of experimental results, and cer- 
tainly does not grow to values of thousands of seconds 
as predicted for a = 0. The validity of our ratchet 
model is quite remarkable in view of its extreme sim- 
plicity; but one must be aware of the following points 
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of disagreement between the experimental and compu- 
tational results (which illustrate the limitations of the 
model): (1) The computed reversal times show weak, 
non-monotonic, dependence on N which is not observed 
experimentally. (2) The largest computed r rev (r rev = 12 
sec for N = 2000) is slightly larger than the experimen- 
tally measured reversal times. (3) The computational re- 
sults for N < 1000 and N > 3000 cannot be directly com- 
pared with experimental results since the corresponding 
reversal times (r rov < 1 sec) fall below the experimental 
resolution. 

The decrease of the computed reversal times for N > 
2400 can be attributed to the "mean field" nature of cal- 
culation of u>3, i.e., to our assumption that (for a given 
N) the detachment of each motor head leads to the same 
energy gain (see Eq. 0]). In reality, the energy change 
upon detachment of a motor depends, in some complex 
manner, on a number of factors such as the positions 
and chemical states of the motors. Motors which release 
higher energy will detach at higher rates, and the de- 
tachment of these "energetic" motors will lead to the re- 
lease of much of the elastic energy stored in the actin 
track. We, therefore, conclude that within the mean 
field approach, the number of disconnecting motors and 
the frequency of detachment events are probably over- 
estimated. This systematic error of the mean field cal- 
culation increases with N, and the result of this is the 
decrease of r rev in this regime, which is not observed ex- 
perimentally. For even larger values of N (N > 5000), 
the model fails because uj^ > loi and the effective at- 
tachment rate outside the gray shaded area in Fig. [7^, 
uj on = u>2 — u>3, becomes negative, i.e., motors detach 
from the track faster than they attach to it. In con- 
trast, for < 1000 < N < 3000, w§ < ui 3 < uj 2 , and 
the effective attachment rate uj on barely changes upon 
changing the model parameter from a = (ref. [2(|) to 
a = 0.0018 (our model). This seemingly minute change 
in Lu on (which, nevertheless, involves a dramatic increase 
in W3) leads to the following non-trivial outcome: On the 
one hand, the fraction of attached motors remains un- 
changed. The data shown in Fig. [7p corresponds to both 
values of a for which the results for N c /N are indistin- 
guishable. On the other hand, the reversal times drop 
by as much as three orders of magnitude (for N = 3000) 
when a is modified from to 0.0018. This spectacular 
decrease in r rev is, therefore, not the result in the change 
in the number of attached motors (since, for each value 
of N, the same fraction of motors is attached for both 
values of a), but rather can be related to the less regular 
manner by which the motors detach from the track. The 
more frequent stochastic detachments of motors from the 
actin track increases the probability per unit time of mo- 



tion reversal. 

5. CONCLUSION 

We have investigated the dynamics of myosin II mo- 
tors on actin tracks composed of small filamentous seg- 
ments with randomly alternating polarities. The absence 
of global polarity leads to a bidirectional relative motion 
between the motors and the tracks. The characteristic re- 
versal time of this motion is of the order of a few seconds 
and exhibits no particular dependence on the number of 
acting motors. Bidirectional motion with macroscopic 
reversal times has been previously observed for NK11 
motors on microtubules and has been attributed to the 
cooperativity of the motors. According to previously pro- 
posed models, the signature of such a bidirectional coop- 
erative motion is the strong exponential dependence of 
T rev on N. The contradiction of this prediction with our 
experimental results can be reconciled by incorporating 
an additional feature into the model, namely, the effect 
of actin stretching by the walking motors. To reduce the 
associated elastic energy, the off rate of motors increases, 
and many of them detach from the track before complet- 
ing a unit step. This effect reduces r rcv considerably and 
eliminates its exponential growth with N. 

Single molecule experiments have led to a dramatic 
increase in our understanding of the structure and dy- 
namics of individual molecular motors. However, many 
biological processes such as muscle contraction, cytokine- 
sis, and the motion of axonemal cilia and flagella, involve 
cooperative action of many motors, which may be af- 
fected by the structure of the underlying track. This con- 
cept is clearly demonstrated in this work dealing with the 
bidirectional motion of myosin II motors on actin tracks 
with randomly alternating local polarities, but without a 
net preferred directionality at the mesoscopic level. This 
unique type of motion is induced by the forces of in- 
dividual motors whose collective effect is manifested in 
macroscopically large reversal times. At the same time, 
the cooperativity of these forces also increases the elastic 
energy of the track, and thereby limits the growth of r rev . 
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